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ABSTRACT 

In the standard fireball model of gamma-ray bursts (GRBs), the fireball starts with an 
optically thick phase. As it expands, the fireball becomes optically thin at some stage. 
The thermal radiation trapped in the originally opaque fireball then leaks out, produc- 
ing a transient event. The appearance of the event is investigated in the framework of a 
homogeneous, spherically symmetric, and freely expanding fireball produced instantly 
by an explosive process without continuous injection of mass and energy. We find that, 
generally, the event has a time-duration shorter than that of the main burst, which is 
presumably produced by the internal shock after the fireball becomes optically thin. 
The event is separated from the main burst by a quiescent time- interval, and is weaker 
than the main burst at least in a high energy band. Hence, the event corresponds to a 
GRB precursor. The precursor event predicted by our model has a smooth and FRED 
(Fast Rise and Exponential Decay) shape lightcurve, and a quasi-thermal spectrum. 
Typically, the characteristic blackbody photon energy is in the X-ray band. However, 
if the distortion of the blackbody spectrum by electron scattering is considered, the 
characteristic photon energy could be boosted to the gamma-ray band. Our model may 
explain a class of observed GRB precursors — those having smooth and FRED-shape 
lightcurves and quasi-thermal spectra. 
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1 INTRODUCTION 

A standard model for gamma-ray bursts (GRBs) has been 
the fireball model (Goodman 1986; Paczynski 1986). In 
this scenario, it is assumed that by whatever a process, 
a radiation-dominated, optically thick, and baryon-poor 
plasma fluid is suddenly produced in a compact volume. 
The radiation field drives the fireball into relativistic expan- 
sion, so that a significant fraction of the initial energy of the 
radiation is converted to the kinetic energy of the fireball 
(Paczynski 1990; Shemi & Piran 1990). At the end of accel- 
eration, the fireball starts to expand freely with a constant 
Lorentz factor > 100 (Kobayashi, Piran & Sari 1999). Later 
on, the fireball becomes optically thin and the thermal radi- 
ation trapped originally in the fireball starts to leak out. As 
the fireball becomes sufficiently large, its kinetic energy is 
converted to the prompt gamma-ray emission through the 
internal shock and the afterglow emission through the ex- 
ternal shock (Rees & Meszaros 1992, 1994; Paczynski & 
Xu 1994). For a review on the fireball model and the in- 
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ternal/external shocks, see Meszaros (2006), Piran (1999, 
2004), and Zhang & Meszaros (2004). 

In many ways a GRB fireball is like a cosmological 
Big Bang (Peebles 1993). Both theories assume that the 
event (the GRB and the Universe) starts with a state of 
a radiation-dominated plasma of very high temperature, in 
which the production of electron-positron pairs is important, 
and the plasma is optically thick. The energy of the radiation 
drives the expansion of the plasma. The plasma and the ra- 
diation cool down as the fireball expands. At some moment, 
the mass density of the plasma becomes low enough so that 
the plasma becomes optically thin to the photons trapped 
in it, and the fireball (of the GRB or of the Universe) un- 
dergoes a transition from an opaque phase to a transparent 
phase. 

However, a GRB fireball also differs from the Big 
Bang. The Universe is homogeneous and isotropic at a very 
high degree. At the time of recombination the fiuctuation 
in the temperature of the cosmic microwave background 
(CMB) radiation relative to the mean temperature is only 
~ 10"^ (Bennett et al. 1996; Smoot et al. 1992). In con- 
trast, the GRB fireball could be highly inhomogeneous and 
anisotropic. Indeed, the fiuctuation in the Lorentz factor of 
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the fireball is required for the internal shock model to work 
(Rees & Meszaros 1994; Paczynski & Xu 1994). The Uni- 
verse is a closed system so that the total energy in it must 
be conserved. However, the GRB fireball interacts with the 
surrounding matter and photons are radiated away from its 
surface. The central engine of the GRB may also continue 
pumping energy into the fireball even after the prompt GRB 
emission has ended (Burrows et al. 2007, and references 
therein). In addition, it is usually assumed that the Uni- 
verse is infinite, but of course the GRB fireball has a finite 
volume. 

A remarkable success of the Big Bang theory has been 
the prediction of the existence of a blackbody CMB of tem- 
perature « 3 K in today's Universe (Gamow 1948a,b; Alpher 
& Herman 1948) and its detection (Penzias & Wilson 1968; 
Smoot ct al. 1992). The CMB is in fact the cooled remnant 
of the primeval fireball — an echo of the Big Bang. Because 
of the expansion of the Universe, the temperature of the 
CMB is redshifted and is hence very low as the CMB pho- 
tons reach an observer of today. Similarly, we expect that 
the GRB fireball has also a remnant of the radiation ini- 
tially trapped in the fireball, and that the remnant would 
have a quasi-thermal spectrum. However, because of the fact 
that the surface of the fireball is assumed to be moving rel- 
ativistically towards an observer who detects the GRB, the 
temperature of the radiation as measured by the observer 
would be significantly boosted by the Doppler effect. 

Because of the fact that the radiation in the fireball 
has a finite energy, tiic remnant event must have a finite 
duration. The event starts when the fireball is still optically 
thick, while the main burst takes place when the fireball 
is already optically thin. Hence, the remnant event must 
occur before the main burst, with a smaller distance from 
the GRB central engine than the main burst. As we will 
see, the remnant event often has a shorter duration than 
the main burst (due to the remnant event's smaller distance 
from the central engine), and is separated from the main 
burst by a quiescent period of time. Since it is produced 
by the emission from the photosphere of the fireball, the 
remnant event should have a spectrum that is dominated by 
a quasi-thermal component. Hence, at least in a high energy 
band, the remnant event should look weaker than the main 
burst. Therefore, the remnant event should be observed as 
a precursor of the GRB. 

The aim of the paper is to quantitatively investigate 
the properties of the GRB precursors as the remnant of the 
radiation initially trapped in the fireball and look for their 
possible observational consequences. 

In at least several cases, precursors of GRBs have been 
unambiguously detected. GRB 030406, a burst that was 
detected out side of the field of view of the International 
Gamma- Ray Astrophysics Laboratory (INTEGRAL), had 
a precursor that occurred ~ 50 s before the main burst 
(Marcinkowski et al. 2006). GRB 041219a, a burst that was 
detected by the INTEGRAL Burst Alert System (IBAS), 
had a precursor that occurred ~ 260 s before the main 
burst (Vestrand et al. 2005; McBreen et al. 2006). GRB 
050820a triggered the Burst Alert Telescope (BAT) onboard 
Swift by a precursor, and Konus/ WIND by the main burst. 
The precursor occurred ~ 200 s before the main burst, and 
the entire duration of the burst is ~ 600 s (Cenko et al. 
2006). GRB 060124, also detected by both BAT/ Swift and 



Komis/ WIND, had a precursor that occurred ~ 500 s be- 
fore the main burst (Romano et al. 2006). The entire du- 
ration is ~ 800 s, making GRB 060124 one of the longest 
bursts. Another interesting case is GRB 061121, detected 
by BAT/ Swift (also by Konus/ WIND and RHESSI), which 
had a precursor that occurred ~ 60 s before the main burst 
(Bellm et al. 2006; Fenimore 2006; Golenetsku et al. 2006; 
Page et al. 2006, 2007). 

These observed gamma-ray precursors have the follow- 
ing characters: (1) The precursor is separated from the main 
burst by a long period of quiescent time; the time separation 
is comparable to the duration of the main burst. (2) The pre- 
cursor is weaker than the main burst, and has a shorter time- 
duration. (3) The precursor often has a smooth and FRED 
(Fast Rise and Exponential Decay) shape lightcurve. (GRB 
050820a is an exception, whose precursor has a lightcurve 
with at legist two peaks.) (4) The precursor has a much 
softer spectrum than the main burst, and in at least two 
cases (GRB 030406 and GRB 041219a) the spectrum can 
be fitted by a blackbody or a blackbody plus a power-law. 

On the other hand, Koshut et al. (1995) have found 
that about 3 percent of the GRBs detected by the Burst 
and Transient Spectrometer Experiment (BATSE) on the 
Compton Gamma- Ray Observatory (GGRO) have precur- 
sors. With a precursor definition different from that of 
Koshut et al. (1995), Lazzati (2005) has found that about 20 
percent of the long-duration BATSE GRBs have evidence of 
precursor emission.^ Soft precursor activities in X-rays have 
been detected in a number of GRBs by Gmga (Murakami et 
al. 1991, 1992) and WATCH onboard GRANAT (Sazonov et 
al. 1998). The precursors detected by Ginga have a thermal 
spectrum with a temperature ~ 1 2 keV. 

On the thermal emission in GRBs, we would also men- 
tion that Ryde (2004, 2005) has found that up to 30 percent 
of long-duration GRBs detected by BATSE/ CGi?0 have a 
spectrum that can be interpreted as combination of a ther- 
mal peak plus a power- law component. But his results refer 
to the prompt gamma-ray emission in the main burst, not 
the precursor emission. 

In this paper, we consider a very simple model for the 
GRB fireball. By assumption, the fireball is homogeneous 
and spherically symmetric, and expands with a constant 
and relativistic speed. The fireball contains a thermal ra- 
diation field, but the energy density of the radiation is not 
large enough to affect the dynamics of the fireball. The rest 
mass and the kinetic energy of the fireball are conserved. 
For the radiation field in the fireball, the change in the 
number of photons is caused only by the emission from the 
photosphere. Hence, the number of photons is an adiabat- 
ically conserved quantity. The last condition requires that 
the electron-positron pairs have already annihilated. Thus, 
in addition, we assume that the opacity in the fireball is 
given by the constant Thompson opacity. 

Our model corresponds to the post-acceleration phase 
of the GRB outfiow. This assumption is justified by the fact 
that the radius of the fireball at the time of photosphere 



^ Lazzati (2005) adopted a definition for a GRB precursor that 
is in favor of weak preeursor emissions and put no limit on the 
time separation between the precursor and the main burst, hence 
he found more precursors than Koshut et al. (1995). 
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emission is much larger than the radius at the end of accel- 
eration (Daigne & Mochkovitch 2002, and Section 5 in the 
present paper). Our model is similar to that of Goodman 
(1986) in several aspects: both models assume a freely ex- 
panding and optically thick fireball in which an amount of 
energy was generated instantly at a beginning time, and the 
observable event is determined by the evolution of the fire- 
ball and the radiation in it. There is no continuous injection 
of energy or mass, unlike the steady wind model of Paczyiiski 
(1986, 1990). However, our model also differs from that of 
Goodman (1986). In the model of Goodman (1986), the fire- 
ball is very hot and the process of pair production is very 
important, and the photosphere was not explicitly included 
in the calculation. While in our model, the fireball is rela- 
tively cool and hence the process of pair production can bo 
ignored (see Section 5 for a justification of this assumption), 
and the photosphere is calculated in details. In addition, 
the model of Goodman (1986) [as well as that of Paczynski 
(1986, 1990)] was aimed to interpret the prompt emission of 
GRBs (the main bursts), while our model is aimed to inter- 
pret the precursor emission of GRBs. This leads to a cooler 
fireball since a GRB precursor usually has a spectrum that 
is softer than the main burst. 

The paper is organized as follows. In Section 2, wc de- 
scribe the geometry and the kinematics of a GRB fireball, 
the kinematics of a radiation field in it, and the structure of 
the photosphere. In Section 3, we describe a formalism for 
calculating the properties of the precursor event of a GRB 
arising from the emission by the photosphere, including the 
luminosity, the blackbody spectrum, and the photon rate 
observed by a remote observer. In Section 4, we present our 
numerical results. In Section 5, wo derive some scaling re- 
lations for the characteristic quantities of the precursor, in- 
cluding the characteristic time scale, the characteristic total 
energy, and the characteristic photon energy. We also give 
brief justification for some key assumptions in our model. 
In Section 6, wo discuss the effect of jet coUimation, the de- 
pendence of our results on the energy band of the detector, 
and the effect of spectrum distortion by electron scatter- 
ing. In Section 7, we summarize the results and draw our 
conclusions. 

In Appendix A, wo present the simplified results for a 
limiting case: a photosphere with a constant expansion ve- 
locity, which applies to the beginning part of the precursor. 

We remark that thermal precursors of GRBs in the in- 
ternal shock model have been previously studied by Daigne 
& Mochkovitch (2002). The thermal emission from a GRB 
photosphere and its effect on the prompt spectrum of a GRB 
have been investigated by Meszaros & Rees (2000); Meszaros 
et al. (2002); Ramirez-Ruiz (2005); Rees & Meszaros (2005); 
Pe'er, Meszaros & Rees (2006); Pe'er et al. (2007); and 
Thompson, Meszaros & Rees (2007). The emission from a 
photosphere of a Poynting flux dominated magnetized fire- 
ball and its relation to the GRB precursor or the prompt 
emission of the main burst has been discussed by Lyutikov & 
Usov (2000); Giannios (2006); and Giannios & Spruit (2007). 
Our model differs from that adopted in their papers. In their 
work, they have assumed a steady wind model for the fire- 
ball, in which mass and energy are continuously injected into 
the wind at the center. And, the expansion of the fireball in 
our model was driven by radiation, not by the Poynting flux 
as in the magnetic models. 



Bianco, Ruffini & Xue (2001) have calculated the ther- 
mal emission from a fireball dominated by electron-positron 
pairs arising from the quantum vacuum polarization process 
around a charged black hole. The results have been applied 
to interpretation of the precursor emission, the main burst, 
as well as the afterglow of GRBs in a unified way (Ruffini 
et al. 2001, 2002, 2005). 

We also remark that thermal precursors produced by 
the shock wave (or jet) breakout of the progenitor star in the 
coUapsar model of GRBs have been studied by MacFadyen 
& Woosley (1999); MacFadyen, Woosley & Heger (2001); 
Ramirez- Ruiz, MacFadyen & Lazzati (2002); and Waxman 
& Meszaros (2003). 

Finally, whenever an observer is referred to in the pa- 
per, we ignore the cosmological effect. That is, we assume 
that the Universe is Euclidean. The correction of the cosmo- 
logical effect is straightforward, which mainly includes three 
factors: cosmological redshift to the observed photon energy, 
time dilation in the observed lightcurvo duration, and that 
the GRB distance appearing in our formulae should be in- 
terpreted as the appropriate cosmological distance. 



2 THE FIREBALL AS A MILNE UNIVERSE 

A freely expanding, homogeneous, and spherically symmet- 
ric fireball is like a Milne universe (see, e.g., Rindler 1977), 
except that a fireball has a finite volume but the Milne uni- 
verse has an infinite volume. 

Let us denote the time in the GRB's rest frame by t, 
and the spherical coordinates in it by {r, 0, </>}. The center 
of the fireball is at r = 0. Then, the Minkowski metric is 

ds^ = -cdt^ + dr^ -f r^dO" , (1) 

where dfi^ = d(P -\- sin^ Qd(^ is the metric on a two- 
dimensional sphere of a unit radius. 
Define the coordinates r) and ^ by 

t = r7Cosh^, r = c77sinh^, (2) 

the metric can then be rewritten as the Milne metric 

ds^ = —c^drf + erf {d^^ + sinh^ f dfi^) . (3) 

The trajcctorjr of a particle with a constant radial ve- 
locity i; = /3c is a straight line defined by 

r = cl3t, (4) 

aside from that 9 = constant and </> = constant. Equation (4) 
demonstrates that, at any moment of constant t, wc have v oc 
r. This profile of expansion velocity is typical for explosive 
events, e.g. supernovae. 

By equation (2) we have that, along the trajectory of 
the particle, 

^ = arctanh/3 , rj = 7~'^t , (5) 

where 7 = (l — (3^^ ^ = cosh^ is the Lorcntz factor of 
the particle. 

Therefore, rj represents the proper time of the particle, 
and ^ measures the spatial velocity of the particle. The spa- 
tial distance in the radial direction on a hypersurface E,, de- 
fined by a constant rj is, by equation (3), / = cr] d^ = crjS, . 

In Fig. 1, we show the spacetime diagram of the fireball. 
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Figure 1. The spacetime diagram of a fireball as a Milne universe. 

Each point represents a two-sphere. The hypersurfacc { = i;ph{v) 
defines the photosphere (Sec. 2.3), within which (the shaded re- 
gion) the photon optical depth is larger than unity. The photo- 
sphere hypersurface becomes spacelike when rj > rjo/ v'2, and ends 
atT] = ?7ph. 



The fireball has an outer boundary at r = R{t) — cPnt, 
where the expansion velocity is vr = cPr, and the Lorentz 
factor isr= 

2.1 The Kinematics of the Fireball 

The comoving spatial volume of a sphere of a coordinate 
radius ^ on the hypersurface E,, is 



Kom(0 = 47rc 



Jo 



sinh^ ^ = TTC^r?^ (sinh 2^ - 2^ - (6) 



For a constant ^, the comoving volume is oc r/' . 

In the non-relativistic limit (3^ ^ <C 1, we have 
Vcom ^ 47rr*/3. In the ultra^-relativistic limit 7 = cosh^ ^ 
1, we have 

V^om « 2n^'cW « —r' • (7) 
7 

The fireball has a uniform comoving mass density p = 
p{r}) oc r)~^. The total rest mass contained in the fireball is 
then 

M = pKom (Cfl) = ^ (sinh 2Cfl - 2^^) , (8) 

where ^r = arccoshF is the value of ^ at the outer boundary 
of the fireball, and the constant reference mass 



Mo = Anpri^c^ 



(9) 



The kinetic energy of the mass contained in a comoving 
volume element dVcom bounded by two neighbored spheres 
of constant radii ^ and ^ + is 

dEK = (7 - l)pcdVcom = (cosh^ - l)Moc sinh^ fd^ , (10) 



which is independent of rj. Hence, the total kinetic energy of 
the fireball is 

En = I dEn 



" / ' 

Ji=o 

= Mod" Q sinh^ sinh 2^^ + ^^fl) 



(11) 



Equation (11) can also be derived as follows. The stress- 
energy tensor of the particle fluid is Tab = pc^driadrit. 
The energy density measured by a rest observer is 
Tab {d/dtY (d/dt)'' = 7^pc^. The total energy contained in a 
hypersurface Ej defined by t = constant (i.e., a three-space 
in the rest frame) is thus 



A /"^ 2 2 2 J Mo 2 - 

= 477 / y pc r dr = / 7 

Jo c io 



3 2, 
r dr , 



(12) 



where the integral is evaluated on St. Substituting r = c/3t 
and rj = 'y~^t into equation (12) and treating f as a constant, 
we get 

E = MoC I -i°l3'dl3 = ^ Moc^ sinh"' , (13) 



= Moc" / 
Jo 



-i''P''dl3= -Mocsmh-" iR 

o 



where j3 = tanh^ and 7 = cosh^ have been used. From 
equations (8) and (13), E — M(? is just the kinetic energy 
i?K in equation (11). 

When /3|e < 1 (P « 1, the non-relativistic limit), we 
have M « MoPr/3 w AirpR^/i and Eyl w Moc^/3|t/10 w 
ZMv^/lQ. As expected, all these results return to the New- 
tonian values. 

When r = cosh^fl ^ 1 {fiR w 1, the ultra-relativistic 
limit), we have 

1 



M « -MoT' , 
and 



Evi 



2 2 



(14) 



(15) 



In the comoving frame the fireball has a uniform den- 
sity. However, in the rest frame of the GRB, the density 
of the fireball increases with radius because of the Lorentz 
contraction. This effect is particularly important near the 
outer boundary of the fireball, where the Lorentz factor is 
dramatically large in the ultra-relativistic case (Fig. 2). 

In the rest frame, the number density of the fireball 
particles is oc 7p oc 7'' on Ef. The total energy density is 
7'^pc^ oc 7"''. Hence, in the rest frame both mass and en- 
ergy are concentrated to a thin spherical shell at the outer 
boundary. If we consider a thin spherical shell defined by 
R — A.R < r < R with a thickness (in the rest frame) 

AR=^«R, (16) 

then by equations (14) and (15), half of the total mass and 
64.6 percent of the total kinetic energy are contained in the 
shell. 



2.2 The Radiation Field in the Fireball 

Assume that the fireball contains a uniform radiation field 
with a comoving energy density Br and a pressure pr = er/3. 
From the conservation of energy, St oc ri~*. The total energy 
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Figure 2. The Lorentz factor of particles in the fireball as 

a function of radius, at any moment of constant t: ■y = 

1^1 — r^/(ct)^j ^ . At the outer boundary of the fireball (the 
dotted line), the Lorentz factor 7 = F = 100. 



of the radiation is not larger than the kinetic energy of the 
fireball, so that the assumption of a freely expanding fireball 
is valid (A < 1, eqs. 24 and 69). 

The total energy of the radiation defined in the comov- 
ing frame is 



Er,coin = erVcom = TTerC 7] (sinh 2^iJ - 2^r) 



(17) 



Since er oc r}~*, we have £r,com oc As the fireball ex- 
pands, the total energy of the radiation decreases, caused 
by the fact that the pressure of the radiation does work. 
However, since the number density of photons rir oc ?7~^, 
the total number of photons in the fireball is conserved.^ 

The stress-energy tensor of the radiation field is Tab = 
(ev + Pv) drjadi^b + (pi/c^) gab, where Qab is the spacetime 
metric. The energy density in the rest frame is then 
Tab (d/dtr {d/dtf = + (7' - 1) Pr = (47^ - 1) e./3. 
The total energy of the radiation measured in the rest frame 
is then 



1 ) CiAnr dr 



(18) 



Using Ct oc T] ^, ri = J ^t, and r = cf3t, we get 
= ^{erv')ch-' I (47^-l)7^/3^d/3. (19) 
cosh^, we can work out the inte- 



Since /3 = tanh^ and 7 ■ 
gral and obtain 

Er = ^ {eiT}*) c^t~^ cosh^fl sinh^ 



(20) 



^ Strictly speaJjing, the number of photons is au adiabatically 
conserved quantity since the emission of photons by the photo- 
sphere reduces the number of photons in the fireball gradually. 




Figure 3. The hypcrsurfacc on which the total energy of ra- 
diation is defined. The energy in the comoving frame, Ei, 

com, IS 

defined on (rj = constant). The energy in the observer's frame, 
Er, is defined on St (t = constant). St and S,, intersect at r = i? 
(the outer boundary of the fireball), so that Er and Er, com are 
related by equation (21). (The diagonal dashed line represents the 
light cone.) 



Hence, Ei cx . 

If we define S,, with the requirement that E,, intersects 
Tit at r = R (i.e., 77 = t/T — f/cosh^ij; see Fig. 3), we can 
compare the Er defined on St and the i?r,com defined on E^. 
We have then the effective Lorentz factor of the radiation 
field 

_ Er 4 sinh^Cfl 



Tr, 



Er, 



3sinh2^H-2^ij 

When r w 1, we have Fr.eff !^ 1, and Er w Er,i 
(47r/3)erJ?^. When F > 1, we have Fr,eff w 2r/3 and 

2, 



:r Er,c 



fer.3 



(21) 



(22) 



where Cr = ei{rj) is evaluated on E^, and Rm ct is the radius 
of the fireball at time t (Fig. 3). 

At any moment of constant t, we have Er\-^^ oc er oc 
Tj^^ cx j^. Hence, 75 percent of the total radiation energy 
(defined in the rest frame of the GRB) is contained in the 
thin shell at the outer boundary of the fireball with a thick- 
ness given by equation (16). 

We can also compare the total energy of radiation on 
St to the total kinetic energy of particles, in the ultra- 
relativistic limit. By equations (7), (8), and (15), we have 

47r 



Ek 



• d3 2 

-pR c 



(23) 



where p is defined on the S^ in Fig. 3. Hence, by equation 

(22), we have 



Ek 



(24) 
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To make the model self-consistent, we must require that 
A < 1. In our model, p oc r?~^ but er oc r;"'*, and so A oc 
r?~^ oc t~^. Hence, if the condition A < 1 is satisfied at some 
moment, it will remain being satisfied afterwards. As rj — » 
we have Cr/pc^ — > oo and A — » oo. Hence, there must exist a 
transition time r/acc defined by A(77acc) ~ 1- When rj < ?7acc, 
the fireball is accelerated by the radiation field. Our model 
only applies to the free-expansion phase of > rjacc . 



2.3 The Photosphere of the Fireball 

To calculate the photon optical depth in the fireball, we need 
to consider the geodesies of a photon. Assume that a photon 
is emitted from a point on a sphere of a comoving coordinate 
radius ^ = at comoving time rji , moving outwards in the 
radial direction. The geodesies of the photon is described by 
r — ri + c{t — ti) in terms of {t,r}, where ti and ri are 
related to r]i and by equation (2). In terms of [r],^}, the 
geodesies can be written as 



7] = fjie^ 



(25 



Assume that the opacity in the fireball, /t, is a constant. 
Then, the optical depth along the trajectory of the photon 
is (Novikov & Thorne 1973) 



(26) 



where dl = crjdi, is radial distance element in the comoving 
frame, and the integral is along the geodesies of the photon. 
Submitting equation (25) into equation (26), we get 



kMo 



[1- 



-2(5h-«i)1 



When ^1=0 and Be ^ 1, we get r 
returning to the Newtonian result. 
When ^1=0 and T > 1, we get 

_ kMo _ kM 



(27) 

(28) 



where Ri = crji is the radius of the fireball at the time when 
the photon is emitted. 

Since r oc rj^'^, at very early time we must have r ^ 
1 and the fireball must be optically thick. We define the 
photosphere of the fireball as a hypersurface determined by 
T = 1. Then, by equation (27), we get the equation for the 
photosphere 



Cph = ^« + ^ In 



1-^ 



(29) 



(30) 



where 

^ //tMoy/' ^ / 3kEk \ 

The photosphere hypersurface is shown in Fig. 1. It 

starts from r; = as a timolike hypersurface, becomes space- 
like when rj > rjo/V^ and ends at = Tjph, where 

1/2 



-■2^R^ 



(31) 



By the time r] = rjph, the radius of the photosphere shrinks 
to zero. At 77 <C 770, we have ^ph « ^r- 
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Figure 4. The fractions /e^ (cq. 35, two upper lines) and /jv, 
(eq. 36, two lower lines) as functions of rj — the comoving time on 
the photosphere. The solid curves are for F = 100. The dotted 
curves are for F = 1000. The maximum of rj is r/ph r;o, defined 
by equation (31) and marked by the dot. At ?7 = J7ph, we have 
= /jv, = 1/8. 



Submitting equation (30) into equation (27), we get 

-2(?„-«i)j _ (32) 



4[1 



In the ultra- relativistic limit F = cosh^j?, ^ 1, we have 
that ?7ph ~ 7?o (1 — l/8r^) « r?o, and r « »?o/^/i when ^1 = 0. 

Although the mass and the energy are concentrated in a 
thin shell at the outer boundary of the fireball (Sections 2.1 
and 2.2), for a photon emitted from the center of the fireball 
the dominant contribution to the total optical depth comes 
from the central region. By equation (27) (with = and 

> 1), the optical depth is t/2 at ^ = In 2^/2 « 0.3466, 
where 7 = 3\/2/4 w 1.0607. This is caused by the fact that 
the mass density decreases quickly as the fireball expands 
and the optical depth is Lorentz invariant. 

For a photon to arrive at the photosphere at a coor- 
dinate radius ^ph and comoving time r?, the photon must 
leave the center of the fireball at time 771 = 776"^^'", 
by equation (25). The photon arrives at a radius ^' at 
time 77' = 7716^ = 776^ -Sph 'p}je total energy of rar 
diation contained in the spheres that have been passed 
by the photon, evaluated along the geodesies of the pho- 
ton, is J eidVcom oc J 77'"^ sinh^ ^'d^' oc J sinh^ ^'d^'. 
The total number of photons contained in the spheres is 

J TlrdKom OC f Sluh^ C d^' ■ 

Let us define 



e = ln 



(l-t-e^p") 



(33) 



on the geodesies of the photon. By equations (29) and (32), 
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the optical depth at ^} is 



(34) 



The fraction of the energy of radiation contained in the 
spheres that have been passed by the photon as it arrives 
at C}i to the total energy of radiation contained in the 
spheres that have been passed by the photon as it arrives at 
the photosphere, is 



/;e-« sinh^^'d^ 
/(f""" e-«' sinh^ f 

3e^ + 6e~^ - e"'^^ - 8 



(35) 



3g€ph _|_ 6e~^pi» - e'^^pi" - 8 ' 
Similarly, the fraction of the number of photons con- 
tained in the spheres that have been passed by the photon 
as it arrives at ^fj, ^ } , to the total number of photons con- 
tained in the spheres that have been passed by the photon 
as it arrives at the photosphere, is 



f^sinh^^'d^' _ sinh2e-2^ 



/q^p'- sinh^ ^'d^' sinh 2Cph - 2Cph 



(36) 



In Fig. 4, we plot fs^ and Jnt as functions of r] on the 
photosphere. The solid curves are for T = 100. The dotted 
curves are for F = 1000. The values of 1 — fs, and 1 — /jv, 
give, respectively, the fraction of the total radiation energy 
and the fraction of the total number of photons contained in 
the spheres that have been passed by the photon on its last 
half journey to the photosphere. When r] < 0.999r;ph, these 
fractions arc respectively about 50 percent and 75 percent. 
As rj rjph, the fractions approach 87.5 percent. 

Inside the photosphere, the radiation has a thermal 
spectrum with a comoving temperature T oc r;^^. The ef- 
fective temperature of the radiation emitted by the photo- 
sphere, measured in the particle's comoving frame at the 
photosphere, can be approximated by 

Teff = Tr"'/" , (37) 

where T is evaluated on the photosphere. Submitting equa- 
tion (34) into equation (37), and letting T — Torjo/ri where 
To is the comoving temperature of the radiation at 77 = »7o , 
we get 

\ -1/2 / 2\ -1/4 

- -1/4 

4e2«K 



(38) 



At ?7 = 77ph (where ^ph = 0), we have t = 1 and Tes = 
To [1 - exp(-25fl)]"^''^ When 77 < r?o, we have r w 'irilh^ 
and Teff « Z-^'^To{r,lrio)-^'\ 



3 THE QUASI-THERMAL PRECURSOR 

In Section 2.3, we have seen that the total optical depth for 
a photon emitted from the center of the fireball is oc r;"^, 



where 77 is the comoving time when the photon is emitted 
(eq. 32). When rj > rjph, where ?7ph is defined by equation 
(31), the whole fireball is transparent to photons. When ij < 
rjph, the fireball has a photosphere, within which the plasma 
is optically thick. The emission from the photosphere has a 
quasi-thermal spectrum. As we will sec, since the emission 
has a finite amount of energy, it produces a transient event 
corresponding to the precursor of a GRB. 

The photosphere is determined by equation (29). From 
this section onward, we assume the ultra-relativistic limit 
F ^ 1. Then, we have ?7ph ~ rjo and e^^ « 2F. Using the co- 
ordinates {t, r}, the photosphere hypersurface is determined 
by the functions of parameter 77 



rph = Fc?7 v'l - 11^ /Vo 



V = Fr? - vVv^o + 



4FVl-'7V%', 



4FVl-W^g/ ' 



(39) 



(40) 



where < 77 < r7ph. 

The velocity of a particle on the photosphere is /3phC = 
fph/iph- The Lorentz factor is 7ph = tph/'?- 

The specific flux density of the radiation emitted by 
the photosphere as measured by a remote observer, whose 
distance to the fireball is much larger than the radius of the 
fireball, is 



-I 



lE„i,^dQohs , 



(41) 



where Bobs is the observed photon energy, Ie^^,^ is the ob- 
served specific intensity of the radiation, and dQobs is the 
element of the solid angle subtended by the image of the 
photosphere on the observer's sky. 

Because of the fact that /bj^^.^ /Sfocai is invariant along 
the path of a photon, where iSiocai is the energy measured by 
any local observer on the path (Misner, Thorne & Wheeler 
1973), we have Ie^^i,^ = Ie^^E'^^^/ E'cm, where Sem is the 
energy of the photon at its point of emission on the pho- 
tosphere as measured by an observer comoving with the 
particle emitting the photon, and /Bem is the specific in- 
tensity measured by that observer. Then, equation (41) can 
be rewritten as 



Fe^ 



I 



g lEetnd^ohs , 



where 



-Eem 7ph (1 - /3ph COS 9) 



(42) 



(43) 



is the Doppler factor of the photon, where 9 is the angle 
between the photon's wave vector and the velocity of the 
particle emitting the photon (Jackson 1999). 

The local specific intensity of the blackbody radiation 
emitted by the photosphere has a Planck spectrum 



Ie^, 



2Ei 



h^C^ exp (Bem/fceTefF) - 1 



(44) 



where fce is the Boltzmann constant, h is the Planck con- 
stant, and Teff is the effective temperature measured in the 
comoving frame (eq. 38). Substituting equation (44) into 
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equation (42), we get 



exp {Eohs/gkBTci 



(45) 



where Ecm ~ -Eobs/g has been used. 

Without loss of generality, we assume that the observer 
is located on the polar axis of the fireball {6 = 0), having 
a distance D from the center of the fireball. At time tph, 
a photon is omitted by a particle on the photosphere. The 
velocity of the particle has an angle d from the polar axis. 
If the photon is emitted in the direction of the observer, it 
will arrive at the observer at time 

iobs = tph - ^ cos 6 , (46) 

where we have shifted the origin of the observer's time so 
that tabs = if the photon is emitted when the fireball has 
a zero radius. To arrive at the observer at the same time, 
photons with different polar angles on the photosphere must 
be emitted at different time. 

In Section 2.3 we have seen that the photosphere hy- 
persurface becomes spacelike when rj > r/o/V^. That is, af- 
ter Tj = tiq/\/2, the photosphere shrinks superluminally. At 
7] = rjo/V^, we have 7ph = tph/n « ^/V^, rph « ctph = 
rcT;o/2. and tph — Vph/c w r/o/2r. As we will see, because 
of the relativistic beaming effect, the dominant contribution 
to the observed radiation comes from a very small region 
on the photosphere around 9 = 0. Hence, by equation (46), 
V = Vo /V2 corresponds to tabs ~ Vo /SF, and 



_ 1 /3kEk\ 



1/2 



(47) 



is a critical time parameter of the model. 

By equations (39) and (40), equation (46) can be trans- 
formed to a cubic equation 



(2m sec' ^) (x' + l) + ( 



a:' -h l) ( 1 + 4r' tan' -]x 



where 



11 = 



0, (48) 



(49) 



Since ^ r} ^ rjph, we have ^ x ^ Xph, where 



(50) 



ph 



Equation (48), which can be solved analytically, has up 
to three real roots of x in the range of ^ a; ^ Xpb, de- 
pending on the values of 9, fj,, and T. From these roots, we 
can calculate the corresponding rj by the inverse of the first 
equation in (49). Hence, we can get the comoving time when 
the photon is emitted from the photosphere 



T?em = Vcm (pi, 9) , 



(51) 



for a given F ^ 1. 

The element of the solid angle, after integration over 
the azimuthal angle 4>, can be written as a function of 9 and 

iobs 



d^^obs — 



27rr 



ph 



D2 



sin 9 cos 9d9 



(52) 



As usual, the azimuthal angle </> spans a range of 0-27r. 
However, the polar angle 9 that contributes to the emissions 



received by the remote observer spans only a range of O-^n 
where (Rees 1966; Bianco et al. 2001) 



• 1 TT 

arccos Pph = arcsm < — 

7ph 2 



(53) 



This is caused by the following fact. Because of the rela- 
tivistic beaming effect, the radiation emitted by a particle 
on the photosphere into the hemisphere of a solid angle 2i: 
around the normal to the photosphere in the particle's co- 
moving frame, spans only a solid angle 27r(l — /3ph) around 
the normal in the rest frame. Photons emitted to the outside 
of this solid angle will hit the optically thick interior of the 
photosphere and be absorbed and hence cannot reach the 
observer (Rees 1966; Bianco et al. 2001). 

Then, substituting equation (52) into equation (45), we 

get 



' ph 



i}{9m-9) sm9 COS 9d9 



exp{Eohs/gkBT^{!) 



(54) 
9) = 1 



where i9(9m — 9) is a step function. That is, i?(( 
if 9 < 9in; = otherwise. 

After solving for r/em as a function of iobs and 9 (eq. 
51), wc can evaluate rph and tph (and hence /3ph and 7ph) 
as functions of tohs and 9 by equations (39) and (40), g by 
equation (43), Teg by equation (38), and 9^ by equation 
(53). Then, the integral in equation (54) can be evaluated. 

We remark that, when equation (48) has nmltiple roots, 
the contribution of each root to the spectrum integral in 
equation (54) (as well as to the luminosity and photon rate 
integrals in eqs. 55 and 56 below) should be summed. 

The luminosity as measured by the observer is then 



L = 4nD^ 



= STTfTsB / g^'ltgr'ph'&iem -9) sin 9 COS 9d9 , (55) 
Jo 

where ctsb = 27r*feB/(15/i'*c') is the Stefan-Boltzmann con- 
stant. It can be checked that when the Tph and T^g are con- 
stants and g = 1, equation (55) gives the standard result for 
the luminosity of a spherical blackbody. 

The lightcurvc of a GRB or its precursor is generally 
not expressed in luminosity. Instead, it is expressed in the 
photon rate, i.e., the number of photons received by the 
observer per unit time. By equation (54), the photon rate, 
which we denote by A^, is 



Es^i^^dEobs 



Tr/2 



N 



47rL> 



/ Fe^^, Jin Eobs 
Jo 

240C(3)crsB 



/2 



g'^TegrlhHOn, - 9) sin 9 cos 9d9 , 



(56) 



where the Riemann zeta function ij(3) w 1.202. 



4 NUMERICAL RESULTS 

At very early time when if /rja <g; 1, the photosphere ex- 
pands with a constant velocity that is close to the speed of 
light, and the formulae for the spectrum, the luminosity, and 
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the photon rate can be simphficd considerably. The results 
in this limiting case are presented in Appendix A. In fact, 
the luminosity and the photon rate can be worked out an- 
alytically. It turns out that the luminosity L is a constant, 
and the photon rate N oc fi^^^ , where the dimensionless time 
H is defined by equation (49). 

Although the spectrum cannot be worked out analyti- 
cally, a remarkable feature of the spectrum is still identified. 
When tobs, EohB, and -Eobs-pEobs scaled respectively by to 
(eq. 47), -Eph.o (eq. A5), and Fq (eq. A8), the spectrum does 
not depend on the value of F, provided that F > 100 (eq. 
A9). As will be seen in this section, this 'similarity' feature 
persists in the late evolution of the spectrum, the luminosity, 
and the photon rate up to a time /i « 10. 

From equations (43) and (46), we have g = Tri/nrjo. 
Then, using the variables x defined in equation (A5) and x 
defined in equation (49) , we can write the spectrum integral 
in equation (54) as 



4 r^' x'fi{x)d{e^,-e) 

Jo exp ij,xx~^/^f2''^{x) 



sin 9 cos 9d9 



where = Eoh^FE^^JFo (eq. A8), 

■^^^^)"(TT^-4fj) ' 

and 



ler^ 



(l + 2^/^/TT^)^ 



(57) 



(58) 



(59) 



Similarly, for the luminosity and the photon rate, we 



have 



and 



31^ 

7n^ 



/ 

Jo 



r/2 



a;*/i (a;)/2"^ (a;)i5(6»i„ - 9) sin 9 cos 9d9 , (60) 



5r^ 



(8v/2- 

7r/2 



■2) 



x''^^fi{x)f^-"''{x)i){9n, - 9)sm9cos9de , (61) 



f-3/4 



where = L/Lo, Lo is the luminosity as /i — > (eq. All); 
and ^fjv s N/No, No is defined by equation (A13). 

Equations (57), (60), and (61) have the forms that are 
convenient for numerical calculations. 

First, we calculate the luminosity as a function of n 
and F. The results are shown in Fig. 5, where the solid line 
is for F = 100, and the dashed line is for F = 1000. At 
small /i, the solid line is almost identical to the dashed line. 
Quantitatively, for /i < 10, the fractional difference in the 
two luminosities is < 10 percent. As /i ^ 0, the luminosity 
L approaches the constant Lo (Appendix A). From /i ~ 0.1, 
the luminosity starts to decay with time. At /i = 1, L drops 
to 0.049Lo- At very large fj,, the luminosity depends on the 
value of F. 

The luminosity drops to zero as approaches F. This is 
caused by the following fact. On the photosphere at any co- 
moving time ?7, the photon emitted by a particle at a larger 9 
reaches the observer at later time. By equation (46), the last 
photon reaches the observer at fobs = tph — rphcos^m/c = 
tph (l — /3ph) = v/lpii! since photons emitted by particles 



■J. -4 



•9- 
< 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



-2-1012 
log M 

Figure 5. Upper panel: the luminosity of the precursor as a func- 
tion of the observer's time. The horizontal axis is /i = iobs/^O- The 
vertical axis is 'i! ^ = L/Lq. The solid line is for F = 100. The 
dashed line is for F = 1000. Lower panel: the fractional difEerence 
in the luminosity: A*/,/*/, = *i(r = 1000)/*/, (F = 100) - 1. 
The dotted line marks the position of A^f/, = 0. 




log ytt 

Figure 6. The fraction of the accumulated energy in the total 
energy (eq. 62). The solid line is for F = 100. The dashed line is 
for r = 1000 (almost indistinguishable from the solid line). 



at ^ > are invisible to the observer. Since 7ph ^ 1, we 
have tobs ^ '7 ^ Vo: hence ^ = toha/to < F. That is, the last 
photon reaches the observer at = F. 

We can define the accumulated energy of the precursor 



© 2006 RAS, MNRAS 000, 1-18 



10 Li-Xin Li 



0.4 - 




12 3 4 



Figure 7. The lightcurve of the precursor. The horizontal axis is 
the observer's time in units of to- The vertical axis is the number 
of photons per unit time, in units of Nq. The solid curve is for 
r = 100. The dashed curve is for F = 1000. 

by 

/■'oba /-M 

E (iobs) = / idtobs = Eo ^Ldn , (62) 
Jo Jo 

where 

Eo = Loto = T-'-Lorjo (63) 

is a characteristic total energy scale. The total energy emit- 
ted by the photosphere is then iStotai = E{tobs = ^to), 
which is equal to 0.426^0 when T = 100, and 0.421i;o when 
r = 1000. 

In Fig. 6 we show the fraction of the accumulated energy 
in the total energy. At n = 1, about 96 percent of the total 
energy has been accumulated. 

The lightcurve defined by ^'at = N/No as a function of 
/i is plotted in Fig. 7, which has a smooth and FRED shape. 
It starts increasing with time as oc fi^^"^, reaches a peak 
at = 0.221, then decreases with time quasi-exponcntially. 
The lightcurve has a long tail of soft photons, which lasts 
until /i = r (sec Fig. 5). These soft photons originate from 
the spacelike part of the photosphere (Fig. 1). 

The width of the lightcurve pulse at iV = iVniax/4 is 
Afi = 0.895. The width at iV = iVmax/20 is A/i = 1.68. The 
maximum of the photon rate varies slightly with F. When 
r = 100, we have A^^ax = 0.363iVo- When F = 1000, we 
have A'max = 0.359A^o- The variation in A^max is at a level of 
1 percent. 

The spectra numerically calculated by equation (57) is 
plotted in Fig. 8, for a sequence of time from = 10""^ to 
H — 10°'^, with the step in log/u being 0.3. As mentioned 
above, for the range of time being of interest (0 < /n < 10), 
the spectrum is insensitive to the variation of F, provided 
that F > 100, -Eobs is scaled by -Eph,o, and tobs is scaled by 
to- In Fig. 8, we plot the spectrum with F = 100 with solid 
lines, and the spectrum with F = 1000 with dashed lines. 




log X 



Figure 8. The blackbody spectrum of the precursor. The hor- 
izontal axis is X = -E'obs/^'ph.Oi which represents the photon 
energy measured by the observer. The vertical axis is = 
-'^obs-'^Boba/-'^0' which represents the energy flux density measured 
by the observer (see the text for details). Different curves corre- 
spond to different observer's time, = tobs/io = 10~^-10°'^ 
(Alog/U = 0.3; from right to left as indicated by the arrow). The 
solid curves are for F = 100. The dashed curves (visually indis- 
tinguishable from the solid curves) are for F = 1000. 

As can be seen from the figure, the two results are almost 
indistinguishable. 

A critical parameter characterizing the energy of the 
photons is the peak spectral energy, iSpeak, which is defined 
to be the photon energy at the maximum of Bobs^Bobs • 
Fig. 9, we plot Speak as a function of time. Since the results 
are not sensitive to the variation of F as we have claimed by 
Figs. 5 8, in Fig. 9 (similarly in the following figures unless 
otherwise stated) we show only the plot for F = 100. We 
see that, as /i — » 0, Speak approaches the asymptotic value 
given by (Appendix A) 

Speak = 4.913 Sph.oM"''" ■ (64) 

Beyond the maximum of the lightcurve (i.e., when /i > 
0.221), the value of Speak decays with time faster than that 
given by equation (64). 

In Fig. 10, we show the spectrum at the maximum of 
the lightcurve (*is,peak = *b at = 0.221, Fig. 7), and 
the spectrum integrated over the time interval of TV/ATmax > 
1/20 

*S,int = / ^EdfX , (65) 
J (11 

where /xi = 0.000345 and /i2 = 1.68. At the low energy end, 
^'s,int oc x^, .just as the instant spectrum at amy moment. 
At the high energy end, *£;,int oc x^^, if we take /ui = 0. 

The integrated spectrum, 'i's/mt, peaks at x = 8.241. 
The spectrum at the maximum of the lightcurve, ^B,peak, 
peaks at X = 10.02. Therefore, the peak spectral energy at 
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Figure 9. The peak energy of the blackbody spectrum as a func- 
tion of time (r = 100). The point marks the peak energy at the 
maximum of the lightcurve (at fj, = 0.221). The dashed line is the 
peaJi energy calculated by equation (64). 



the maximum of the hghtcurve is 

-^^peak.max 10.02 £;ph,o. (66) 
The peak spectral energy of the integrated spectrum is 

■E'peak.int = 8.241 Sph,o ■ (67) 

5 CHARACTERISTIC QUANTITIES 

Our model has three characteristic quantities: i?ph,o (eq. 
A5), a photon energy scale charactering the peak energy 
of the spectrum; to (eq. 47) , a time scale characterizing the 
duration and the 'pulse width' of the lightcurve; and Eq (eq. 
63) , an energy scale characterizing the total energy emitted 
by the photosphere. 

The ratio Eg / Ek determines the dynamical importance 
of the radiation field to the fireball. As mentioned in Section 
2.2, the importance of the radiation field for the dynamics 
of the fireball is determined by the parameter A defined by 
equation (24), where Cr is the comoving energy density of 
the radiation, and p is the comoving mass density of the 
fireball. Since fir — aT*^ = aTQ{ri/rio)^'^, where T is the co- 
moving temperature, and a — 4(tsb/c is the radiation den- 
sity constant, by equations (7), (8), and (15) we have 



X: 



3Sk V Vo 



(68) 



By equations (A6), (All), and (63), we have then 



6Lo?7o 




6eo 




7r£;K 




7 





where 



Eq 

Ek 



(69) 



(70) 




log X 

Figure 10. The integrated spectrum of the precursor ('I'B,int, 
integrated over the time interval of N/Nmux > 1/20, eq. 65; 
solid curve), and the spectrum at the maximum of the lightcurve 
{'^ E , peaky at // = 0.221; dashed curve). The dotted line is the 
asymptotic spectrum at high energy end when the spectrum is 
integrated from /i = 0, ^'B,int (^q. A16). 



Hence, a,t rj — rjo, we have A — 6eo/7. For the model of 
a freely expanding fireball to be a good approximation, we 
must require that eo <C 1. 

From equations (63), (A5), (All), and the definitions 
of to (eq. 47) and .Ro (eq. A6) , we can derive that 



iph,0 



3/4 



-1/2 



(71) 



If Eo, Eph,o, and to can be measured, then equation (71) can 
be used to estimate the value of the Lorentz factor P. 

Since r > 1, equation (71) leads to a constraint on Eq, 

to, and -Eph.o 



-3/4 



(72) 



As we have assumed, the fireball does not contain 

electron-positron pairs and the opacity in the fireball is dom- 
inated by the Thompson electron scattering opacity. Then 
we have k = 0.20(1 -|- X) cm^ g~^, where X is the mass 
fraction of hydrogen. It is likely that long-duration GRBs 
arise from the core-collapse of massive stars which have lost 
their hydrogen envelopes (Woosley & Bloom 2006), hence it 
is reasonable to assume that X = 0. Then we have k = 0.2 
cm^ g"^ 
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to 



Then, by equation (47), we have 

, 1/2 

1.719 s ' 



0.2cm2g-i J \ 1054eri 



Ek 



1/2 



p s -5/2 



By the definition of eo, equation (70), we have 

Ek 



Eo = W^^eri 



Vo.oi 



1054erg 



(73) 



(74) 



By equations (47), (70), and (71), we have 

1/4 / \ -3/8 



Sph,o « 1.381 kev(^ 



0.2cni2g-i 

/ £k y'^V r yv« 

^ \^1054erg J \WqJ 
By equations (75) and (A5), we have 

1/4 



(75) 



/cbTo « 0.018 keV 



f— V 

Vo.ou 



0.2cm2g- 



-3/8 



Ek \-'''/ r x3/« 



10-'4cre 



100 ) 



(76) 



The very small value of ksTo justifies our assumption that 
the fireball does not contain pairs. ^ 

By equations (47) and ( A6) , the characteristic radius of 
the fireball when the thermal precursor emission takes place 
is 



Ro = r^cto w 5 X 10"cm 



1/2 



0.2cm2g-i 

1/2 . T. . -1/2 



(77) 



io54ergy viooy 

The radius of the fireball at the end of acceleration, defined 
by A = 1 (t? = Tjacc), is 

4>^^°"^-((St) (o.2em2g-i'"" 

1/2 , , _l/2 

(78) 



Rax:c — -y--Ro 



^ \^1054erg J \Wo) 



For the model of a freely expanding fireball to be self- 
consistent, we must require that Jiacc <C Ro, i-e. eo ^ 1- 

The peak spectral energy is related to i5ph,o by equa- 
tions (64), (66), and (67). The duration of the precursor (de- 
fined by the width of the lightcurve pulse aX N = Nmax/^O) 
is w 1.7to- The total energy emitted in the precursor, is 
w 0.42So. 

Generally, the characteristic photon energy predicted by 
our simple model with the assumption of a blackbody pho- 
tosphere emission is too small for explaining the precursors 
detected by the gammar-ray detectors, but is in agreement 



3 Even if -Eph.o lias a value as high as ~ 260 keV which is likely 
the case when the blackbody spectrum is distorted by electron 
scattering (Section 6.3), we have fceTo ^ 3 keV and the effect 
of pair production is also unimportant at least during the stage 
when the majority of the precursor photons are emitted. 



with the precursors detected by the X-ray detectors (see In- 
troduction; for further discussions see Sec. 6.3). 



6 DISCUSSIONS 

6.1 The Effect of Jet Collimation 

In our model we have assumed that the fireball is a sphere 
with a total solid angle of Att. In reality, GRBs are likely to 
be coUimated. That is, the fireball may have a cone shape 
and spans a solid angle f2 < 47r. Now we check how this 'jet 
collimation' effect affects our results. 

Assume that the fireball spans a solid angle O = ^ttlo, 
a; < 1. The observer is on the polar axis 9 = 0, and the polar 
boundary of the fireball (which is now a cone) is at = 0c < 
7r/2, where 6c = arccos(l — 2a;). In addition, we assume that 
6c 3> T"^, so that the effect of 'jet collimation' starts to 
affect the shape of the lightcurve only at very late time when 
the Lorentz factor of the photosphere (7ph) drops to a value 
~ 1/dc, just as the observation of the 'jet break' signature 
in GRB afterglow lightcurves (Sari, Piran & Halpern 1999; 
Friedman & Bloom 2005, and references therein). Hence, we 
have < o) < 1. 

With jet collimation, the comoving volume of the fire- 
ball in equation (7) should be replaced by 

Vcom ~ 2TTUjT^d\^ OC LO , (79) 

where we have let 7 = F. The kinetic energy Ek is still 
related to the total mass M by equation (15), but equation 
(14) should be replaced by 

M « |MoF' , (80) 

where AIq is defined by equation (9). 

The optical depth is still given by equations (27) and 
(32), and 770 is still defined by equation (30). Then, for the 
characteristic time in equation (47), we should replace Ek 
by Ek/u>. This is equivalent to say that equation (47) still 
holds if Ek is interpreted as the isotropic-equivalent kinetic 
energy. 

The photosphere is still defined by equation (29), or 

equivalently, by equations (39) and (40). The effective tem- 
perature of the photosphere is still given by equation (38). 

Equation (48) is unchanged, except that is only al- 
lowed to vary in the range of 4 < 0c- 

The specific flux density measured by the observer is 
still given by equation (54), but the following replacement 
to the integral over 6 should be made 



Jo Jo 



de . 



(81) 



The luminosity defined by equation (55), after the re- 
placement in equation (81), should be interpreted as the 
isotropic-equivalent luminosity. The true luminosity of the 
photosphere should be cjL. The same conclusion holds for 
the photon rate in equation (56). 

Because of the relativistic beaming effect, when 7ph » 1 
the contribution of the photosphere to the observed luminos- 
ity, the photon rate, and the spectrum comes from a small 
region of < 7pj^^ <g; 6c on the photosphere. Hence, the effect 
of the boundary of the jet becomes important only after 7ph 
drops to a value smaller than 7ph,br ^ Oc^ . For Oc » F"^, 
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Figure 11. The flux-averaged polar angle on the photosphere 
(solid curves), defined by equation (82). The curve ends at /i = F, 
since after /x = T no photons emitted by the photosphere can 
reach the observer. The dashed lines show the maximum polar 
angle of the conc-shapc fireball, Be- The model of the spherical 
fireball applies only when the value of (6) is below the dashed line 
(i.e., when {0) < 9c). If the photosphere is a static and uniform 
sphere, we would have {6) = 7r/4. 

we have 7ph,br ^ T, and hence the time n = /Xbr when 7ph 
drops to 7ph,br would be S> 1. 

Wc can define a flux-averaged polar angle on the pho- 
tosphere by 

(^0) = / ^ 6l/T5f rpht? (61m - 6l)sin6lcos6ld6l , (82) 

^ Jo 

for the model of the spherical fireball. If the photosphere is 
a static and uniform sphere (i.e., g = 1, Tea and Vph are con- 
stants), we would have (d) = 7r/4. If the contribution to the 
luminosity comes from a small region of ^ <g; 1 on the pho- 
tosphere (which is the case when the photosphere expands 
relativistically), we would have (^) <C 1. If the contribution 
to the luminosity comes from a narrow ring region around 
6 = 6q on the photosphere, we would expect (9) ~ ^o- 

In Fig. 11, we show the value of {9} as a function of 
time (solid curves). When fj, ^ 1, the photosphere has a 
Lorentz factor 7ph Ri F 1, the radiation received by the 
observer comes dominantly from a region of < on the 
photosphere, hence we have {9) ~ F"^ <C 1. As /it — > /Xmax = 
F, the radiation received by the observer comes dominantly 
from a region of n/4 < 6 < 7r/2. 

The dashed lines show the value of ^c, the maximum 
polar angle of the cone-shape fireball. When (9) > 6c, the 
observer starts to see the boundary of the cone and the 
spherical model does not apply any more. As a result, the 
lightcurve starts to decay faster with time than that pre- 
dicted by the spherical model. The transition time /itr (or, 
the 'jet break' time, as commonly called for the GRB af- 
terglow) is approximately determined by (9) = 9c, i.e. the 
intersection of the solid curve and the dashed line. When 
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Figure 12. The lightcurves detected by detectors with different 
lower limit on the photon energy. Top to bottom: Bmin/Sph^o = 0, 
1, 5, 10, 15, and 20. 

9c S> F~^ as we have assumed, the transition time /itr is 
much larger than the duration of the event. 

Particularly, when -C 1, the effect of jet coUimation is 
negligible. The results in Appendix A are unchanged. i?ph,o 
is still given by equation (A5), Lo is still given by equation 
(All), and A^o is still given by equation (A13). However, as 
explained above, Lo and No, as well as the Ek in the defi- 
nition of rjo (eq. 30) , should be interpreted as the isotropic- 
equivalent quantities. 

Then, i5ph,o is still the characteristic photon energy, 
while Eo = Loto is the isotropic-equivalent characteristic 
total energy. Equations (70)-(77) are then unchanged, if the 
Eyl and Eo are interpreted as the isotropic-equivalent kinetic 
energy and the isotropic-equivalent characteristic total en- 
ergy, respectively. 



6.2 Dependence of the Lightcurve and Spectrum 
on the Energy Band of a Detector 

Every detector has a finite range of energy that the detector 
is sensitive to. For example, the BAT on Swift covers an en- 
ergy range of 15-150 keV, while the BATSE on CGRO cov- 
ers an energy range of 20-1900 keV. The observed lightcurve 
and the spectrum of a GRB precursor depends crucially on 
the range of energy of the detector. 

We consider only the effect of the lower limit on the pho- 
ton energy. For a quasi-thermal spectrum the flux density 
decays quickly at the high energy end, the upper limit on 
the photon energy is not important provided that it exceeds 
the peak spectral energy at the maximum of the lightcurve.* 
Hence, we take the lowest energy to be i^min, and allow the 

■* This could be a problem for Swift since the BAT on it has a 
very small maximum energy. 
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Figure 13. The observed peak photon energy of the spectrum at 
the maximum of the lightcurvc (upper panel) and the width of 
the observed Ughtcurve pulse (at N = Afmax/SO; lower panel) as 
functions of the lower energy limit of the detector. The dashed line 
in the upper panel shows the lower energy limit of the detector. 



highest energy to be infinity. Then, the observed photon rate 
is 

^4 



N 



N_ _ _ 

iVo ^ 36 



(83) 



where Xmin = -BminZ-E'ph.O- 

The hghtcurves calculated for Xmin = 0, 1, 5, 10, 15, and 
20 are shown in Fig. 12. As Emin increases, the maximum 
of the lightcurve shifts to earlier time, and the width of 
the lightcurvc pulse decreases. This is caused by the fact 
that the energy of the photons emitted by the photosphere 
decreases with time (Figs. 8 and 9). The maximum photon 
rate decreases as Emin increases. 

In Fig. 13, we show the observed peak photon energy of 
the spectrum at the maximum of the lightcurvc (estimated 
by eq. 64) and the width of the observed lightcurve pulse at 
N = iVmax/20 as functions of Emin- The pealc photon energy 
increases with Emin, while the width of the Ughtcurve pulse 
decreases with iSmin- 



blackbody with temperature ~ 106 ± 13 keV, corresponding 
to a peak photon energy « 416 ± 51 keV (in the observer's 
frame; Marcinkowski et al. 2006). However, the pealc photon 
energy of the integrated spectrum predicted by our model 
is typically ^ 8i5ph,() ~ 10 keV, which appears to be con- 
sistent with that of the soft X-ray precursors detected by 
Ginga (Murakami et al. 1991, 1992) and WATCH/ Gii^iV^T 
(Sazonov et al. 1998). 

In our model we have assumed that the emission from 
the photosphere is a perfect blackbody. However, when elec- 
tron scattering dominates the opacity, thermal ccjuilibrium 
exists only when TgTTh > 1, where Tg is the free-free opti- 
cal depth, and vrh is the Thompson optical depth (Felten & 
Rees 1972; Shakura 1972; Shakura & Sunyaev 1973; Novikov 
& Thorne 1973; Sunyaev & Shakura 1974). At the surface 
of the fireball, the spectrum of the outgoing radiation is dis- 
torted. The energy flux measured in the comoving frame is 
(Novikov & Thorne 1973) 



Q = 1.54 X 10"*ne''^: 



n9/4 

eff 



(84) 



in cgs units, where Ue is the number density of electrons, and 
Teff is the effective temperature. Equation (84) is in contrgist 
to the standard formula Q = asBT^g for the blackbody. 

The characteristic radius for the emission from the pho- 
tosphere is the Ro given by equation (77). Let us define the 
effective temperature of the standard blackbody spectrum 
in the comoving frame, Teff,o, by Lq = 4a:asB.T^gfiRo^^ , 
where Lq is related to the characteristic total energy Eq by 
equation (63). Similarly, we can define the effective tem- 
perature of the distorted blackbody spectrum, Teff,!, by 
Lo = 47rQi^o^^ where Q is given by equation (84). Then, 
we find that 



23.5 
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(85) 



where we have assumed that k = 0.2 cm^ g ^ and Ue = 

p/2mH, where mn is the atomic mass unit. 

If we take -Bpeak ~ 8-Bph,o, then equations (75) and (85) 
imply that 



^peak 



\ s 4/9 

259 keV ( — 1 



(a 



01 y 



lOB^erg 



-1/9 



(loo) 



4/3 



(86) 



6.3 Distortion of the Blackbody Spectrum by 
Electron Scattering 

Even if the effects of jet coUimation and the finite en- 
ergy band of the detector are considered, the characteris- 
tic photon energy predicted by the model is still signifi- 
cantly smaller than that of the GRB precursors detected 
by gamma-ray detectors. For example, the spectrum of the 
precursor of GRB 041219a can be fitted by a blackbody plus 
a power law. The temperature of the blackbody is 46 ± 9 
keV, corresponding to a peak photon energy ~ 180 ±35 keV 
(in the observer's frame; McBreen et al. 2006). The spec- 
trum of the precursor of GRB 030406 can be fitted by a 



for the distorted blackbody spectrum. The peak spectral 
energy given by equation (86), which is defined in the ob- 
server's frame, appears to be in agreement with that of the 
precursors detected by gamma-ray detectors. 

Another effect that may help shifting up the photon 
energy is that the blackbody emission from the photosphere 
may be nonthermalized by some scattering process and have 
a nonthermal spectrum. For example, Thompson (1994) has 
shown that the scattering of the adiabatically cooled thermal 
photons by the magnetohydrodynamic (MHD) turbulence 
in the photosphere can boost the thermal photons up to a 
larger peak energy and lead to a spectrum that have a broken 
power-law or a Band function (Band et al. 1993) shape. 
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7 SUMMARY AND CONCLUSIONS 

Wo have studied a simple model for the precursor event of a 
GRB, where the precursor is produced by the emission from 
the early photosphere of a fireball. The fireball is assumed 
to be homogeneous, spherical, and expand freely and relar 
tivistically. The fireball contains a constant amount of rest 
mass, and has a constant total kinetic energy. It also con- 
tains a radiation field, which drives the acceleration of the 
fireball expansion at the very early stage, but simply cools 
down with the expansion at the late free expansion stage. 
These assumptions make our model different from that has 
been adopted in the literature for studying the photosphere 
emission of a GRB, which usually assumes a steady wind 
with mass and energy injection at the center (Daigne & 
Mochkovitch 2002; Mcszaros & Rees 2000; Mcszaros et al. 
2002; Ramirez-Ruiz 2005; Rees & Mcszaros 2005; Pe'er et 
al. 2006, 2007; Thompson et al. 2007). 

The dynamics of the fireball is the same as that of a 
Milne universe, except that the fireball has a finite size and 
contains a finite amount of mass and energy. In the same way 
as the Universe, at very early time the fireball is opaque to 
the photons contained in it, but later it becomes transparent 
(Fig. 1). The photosphere is defined to be the hypersurface 
where the total optical depth is unity. If the fireball had an 
infinitely large radius, then as in the Milne universe, the pho- 
tosphere were just a spacelike hyperbola. However, because 
of the fact that the fireball has a finite radius, the hypersur- 
face of the photosphere is timelike when < rio /\/2. It is the 
emission from this timelike part that makes the dominant 
contribution to the precursor emission detected by a remote 
observer, since this timelike surface moves towards the ob- 
server with a relativistic speed. However, the spacelike part 
of the photosphere produces a long tail of soft emission. 

The luminosity, the photon rate (the lightcurve), and 
the blackbody spectrum of the precursor are calculated. The 
results axe presented in Figs. 5-10. A remarkable feature of 
the results is that, when the quantities are properly scaled 
and the time is not too large, the results are almost invariant 
with respect to the Lorentz factor of the fireball. That is, 
if the observer time is in units of to = ??o/r (eq. 47), the 
luminosity is in units of Lq (eq. All), the photon rate is 
in units of Nq (eq. A13), and the spectrum -BobsfEots ™ 
units of Fq (eq. A8), then the luminosity, the photon rate, 
and the spectrum are insensitive to the variation of F for 
M = toha/to ^ 10, provided that F > 100. At very late time, 
the results depend on the value of F. 

Thus, a critical time scale for the precursor is to (cqs. 47 
and 73), which is ~ 2 s for typical parameters. The luminos- 
ity starts with a constant L = Lo, and drops to L 0.051/0 
at tobs = to- The total energy emitted by the photosphere 
is ~ 0.42L()t(). About 96 percent of the total emitted energy 
has been accumulated at tobs = to- 

The lightcurve (defined by the photon rate) has a 
smooth and FRED shape (Fig. 7). It begins as A'^ oc t^(^^, 
peaks at tobs ~ 0.22to, and then decays quasi-exponentially. 
The width of the lightcurve pulse at A" = Amax/4 (Amax is 
the maximum photon rate) is Atobs ~ 0.9to. The width at 

N = Amax/20 is Atobs w 1.7to. 

The blackbody spectrum at different moment of ob- 
servation is shown in Fig. 8. The peak spectral energy 
£^peak, defined to be the photon energy at the maximum of 



B^obs-FB^bs > decays with time (Fig. 9). Before the maximum 
of the hghtcurve (tobs < 0.22to), the peak spectral energy 
evolves according to, approximately, Speak oc t~^/'^ (eq. 64). 
After the maximum, the peak spectral energy decays with 
time with a rate faster than that given by equation (64) 
(Fig. 9). At the maximum of the lightcurve (tobs ~ 0.22to), 
the spectrum peaks at S'peak = 10.0Sph,o, where Ep^fi is the 
characteristic photon energy defined by equation (A5). The 
spectrum integrated over the time interval of A > Ainax/20 
peaks at Speak = 8.24£;ph,o (Fig. 10). 

The characteristic quantities describing the precursor, 
which include the characteristic time, the characteristic to- 
tal energy, and the characteristic photon energy are sum- 
marized in equations (73)-(75). For typical parameters, the 
characteristic time is shorter than the duration of a typical 
long GRB (the main burst), and the characteristic photon 
energy is in the X-ray band. Because of the thermal nature 
of the spectrum, at least in the high energy band the to- 
tal energy is expected to be smaller than the main burst. A 
constraint on these three characteristic quantities is given 
by equation (72), which does not depend on the value of the 
opacity in the fireball. 

The characteristic radius of the thermal emission, Rq, 
is given by equation (77), which is 10^"'-10^^ cm for typical 
parameters. The internal shock must not occur until the 
fireball becomes optically thin, since otherwise it would be 
difficult to explain the nonthermal spectrum of the GRB. 
Hence, Rg — the radius where the internal shock occurs — 
must be > _Ro. The time separation between the precursor 
and the main burst is then tsep ~ Rs/2r^c. 

The main burst produced by the internal shock has 
a duration T ~ A/c, where A is the width (measured in 
the rest frame of the GRB) of the outflow that contributes 
to the prompt gamma-ray emission (Piran 1999). In our 
model, the fireball has an expansion velocity v <x. r, and 
the majority of mass and energy is confined in a thin layer 
at the outer boundary of the fireball (Sec. 2.1). The thick- 
ness of the layer — which contains 50 percent of total mass 
and 65 percent of total kinetic energy — is given by equa- 
tion (16). Letting R = Rs, we have then A w Rg/2r^ and 
T fs Rs/2V^c ~ tscp. This relation between the duration of 
the main burst and the time separation between the main 
burst and the precursor appears to be consistent with the 
observation (see Introduction). 

We have assumed a spherical fireball. If the GRB out- 
flow is collimatcd and the resulted fireball has a conflgu- 
ration of a cone spanning a solid angle Anuj that satisfies 
<t. oJ < 1, the results are still valid until a time is 
reached when the flux-averaged polar angle of the photo- 
sphere increases to a value ^ Oc — arccos(l — 2uj) (Fig. 11). 
However, in the results, L and N should be interpreted as the 
isotropic-equivalent luminosity and the isotropic-equivalent 
photon rate respectively. Equations (70)-(77) still hold, but 
Ek and Eq should be interpreted as the isotropic-equivalent 
kinetic energy and the isotropic-equivalent characteristic to- 
tal energy of radiation, respectively. This fact is caused by 
the relativistic beaming effect, which implies that when the 
photosphere expands relativistically the dominant contribu- 
tion to the spectrum and the luminosity observed by a re- 
mote observer comes from a small region on the photosphere 
that spans a solid angle ~ F~^. 

The dependence of the observed lightcurve and spec- 
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trum on the finite energy range of a detector lias also been 
discussed (Sec. 6.2). Overall, as the lower limit of the de- 
tector energy increases, the duration of the lightcurve de- 
creases and the maximum of the lightcurve shifts to earlier 
time (Figs. 12 and 13). The peak photon energy of the spec- 
trum increases with the lower energy bound of the detector 
(Fig. 13). 

The photon energy predicted by our simple model seems 

to be consistent with the X-ray precursors of GRBs detected 
by G'lnga and WATCU/GRANAT (Murakami et al. 1991, 
1992; Sazonov et al. 1998), but is too small for the precur- 
sors detected by gammarray detectors (Cenko et al. 2006; 
Marcinkowski et al. 2006; McBreen et al. 2006; Romano et 
al. 2006; Page et al. 2007). However, in our model wo have 
assumed that the emission from the photosphere has a spec- 
trum that is perfectly blackbody. When electron scattering 
dominates the opacity, the blackbody spectrum is distorted 
(Felten & Rees 1972; Shakura 1972; Shakura & Sunyaev 
1973; Novikov & Thornc 1973; Sunyaev & Shakura 1974). 
We have estimated this effect by adopting a modified fiux 
equation in equation (84) and found that the distorted spec- 
trum could have a pealc spectral energy of ~ 260 keV, which 
is consistent with the hard spectrum of the observed gamma- 
ray precursors (Sec. 6.3). Scattering of the soft blackbody 
photons by e.g. the magnetic turbulence in the fireball may 
also lead to a nonthermal spectrum (Thompson 1994). A 
detailed study on these issues is left for our future work. 

A strong assumption in our model is that there is no 
continuous injection of energy and mass, and hence the to- 
tal mass and the kinetic energy of the fireball are conserved. 
A big benefit from this simplified assumption is that the 
shape of the lightcurve and the duration of the omission can 
be predicted, in contrast to the steady wind model where 
the lightcurve and the duration depend on the rate of en- 
ergy and mass injection as a function of time. On the physics 
side, this type of models with 'instant release of energy' may 
be relevant to the precursors that have smooth and FRED 
shape lightcurves. More complicated models with continu- 
ous injection of mass and energy may be relevant to the 
precursors that have complex lightcurves. 

Although the X-Ray Telescope (XRT) on Swift has dis- 
covered that flares arc common in early X-ray afterglows of 
GRBs and it seems that the most likely cause of the flares is 
the late-time activity of the GRB central engine (Burrows 
et al. 2007, and references therein), the time separation be- 
tween two adjacent flares are usually much longer than the 
duration of the precursor event predicted by our model and 
hence the assumption of no continuous energy and mass in- 
jection remains valid. Of course the fireballs or expanding 
shells related to the late-time activity would also have photo- 
sphere omissions. However, the late-time photosphere omis- 
sions would overlap with the prompt main GRB emissions 
or the afterglows (Meszaros & Rees 2000; Meszaros et al. 
2002; Ramirez-Ruiz 2005; Rees & Meszaros 2005; Pe'er et 
al. 2006; Thompson et al. 2007) and hence are hard to detect 
[see, however, Ryde (2004, 2005); and Pe'er et al. (2007)]. 

In summary, our model implies the existence of a quasi- 
thermal precursor of a GRB. The precursor is a remnant of 
the thermal radiation contained in the fireball during its ini- 
tial acceleration phase, just like that the CMB is a remnant 
of the radiation in the Big Bang. The precursor emission 
has a quasi-thermal spectrum, and a smooth and FRED 



shape lightcurve with a duration ~ 1-5 s. If the distortion 
of the blackbody spectrum by the electron scattering pro- 
cess is not important, the radiation observed by a remote 
observer is in the X-ray band with a peak photon energy 
~ 1-10 keV. If the distortion by electron scattering is im- 
portant, the peak photon energy could be of several hundred 
keV and hence be in the gamma-ray band. Under some con- 
ditions the soft thermal photons may also be nonthermalized 
and have a nonthermal spectrum. Although in reality the 
situation may be much more complicated, our simple model 
may provide a reasonable interpretation for at least a class 
of GRB precursors — those having smooth and FRED-shape 
lightcurves and quasi-thermal spectra. 
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APPENDIX A: PHOTOSPHERE WITH A 
CONSTANT EXPANSION VELOCITY 

When rf /r^Q <ti 1, the photosphere defined by equation (29) 
has a radius ^ph « £,r, constant 7ph = cosh^ph ~ T » 
1, and constant dm ~ l/F ^ 1. That is, the photosphere 
expands with a constant and ultra-relativistic velocity. In 
this Appendix, we derive the spectrum, the luminosity, and 
the photon rate of a photosphere in this limiting case. 

When rj^/riQ 1, we have tph ~ Trj, rph w ii w Fa?, 
and /3ph w /3ij w 1 — l/2r^. Then, by equation (43), we have 



9 



1 



(Al) 



r(l- /3phCos6») 
By equations (43) and (46), we have tobs = v/d and hence 

■no r 

Then, by equations (38) and (A2), we have 

1/2 



(A2) 



-1/4 



J-eB ■ 

and 

-Eobs ^ 

gkBTeff 
where 

-Eobs 



Vo 



1/2 



X 



(A3) 



(A4) 



(A5) 



-Bph.O 

Let us define 

Ro = Tcr^o ■ (A6) 
Then, Vph/Ro ~ v/Vo ~ m/^- equation (54) we have 



WasBToRo 2 2 4 

-r M X 



p-1 



g sm tl cos t 



Define s = g/T and 



Fo 



Fo = 



47r£)2 ' 



(A7) 



(A8) 



where Lq is the constant luminosity in the limit rj^ /t/q <C 1 
given by equation (All) below. By equation (Al), we have 
cos6> « /3ph [l - l/(r^s)] and dcos6' « p-^^ds/iV^s^). 
When ^ = 0, we have s = 2. When 6 ^ l/F, we have s w 1. 
Hence, equation (A7) can be rewritten as 



45 2 4 



ds 



exp[(,./s)V2x]-l 



(A9) 
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In equation (A9), x s-^d /i appear in the combination 
lj}^^X- Hence, ^'b is invariant under the transformation 

/ \ 1/2 

M-m', X^x'=xiA] ■ (AlO) 

That is, the value of '^e at energy x and at time /u, is equal 

to the value of 'I'e at energy fi^^^x. and at time // = 1. When 
H = 1, ^'b(x) peaks at x = 4.913. Hence, at riw /i, 4'b(x) 
peaks at Xpeak = 4.913 

Similarly, for the luminosity in equation (55), we have 

L w ^crsBToi?o / 3^ sin ^ cos 
^ Jo 

« Lo = ^asBT^Rlv\ (All) 
9 

which is a constant. 

For the photon rate in equation (56), we have 



_ 80-3^/^C(3)c7SB ^3 ^2 
^ ~ =3^ 



p-1 

X / 5 



r-1 / \ 1/2 

,3 / 



= ATom'/', (A12) 
where 

Equations (All), (A9), and (A12) approximate the nu- 
merical results presented in Figs. 7-9 very well for /i < 0.1, 
with fractional errors < 5 percent. 

The integration of ^Pb over ji is 



45 _2 w^dw ,,,,, 
^^Ed^l = —jx ds = , (A14) 



exp Ww/s — 1 



where w = nx^ . 

As X — *■ 0, we have 

jT" ^,d„ « (2^/^ - 1) M^/^x^ oc x^ . (A15) 

As X — *■ we have 



onr»— 2 



*,d,.«^X-Vx-^ (A16) 
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